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Abstract 

The mathematical models of many 
dynamic systems of interest in the aerospace 
industry are inherently complex and of high 
order. Rather than grapple with the full- 
complexity models for such systems, the control 
designers often elect to derive low-order, reduced- 
complexity models for which a control system 
can be designed. The subject of model-order 
reduction has received a significant amount of 
attention in control engineering literature. This 
paper describes and demonstrates a novel active- 
control methodology that can be used to design a 
control system that automatically forces the 
original system to behave like a chosen reduced- 
complexity model by treating the effects of the 
original complexity-related terms as disturbances 
acting on the reduced-complexity system. The 
method is demonstrated by several worked 
examples involving both linear and nonlinear 
systems, supported by simulation results. 

Introduction 

The difficulty in developing high-order 
models for complex systems in the aerospace 
industry has been facilitated by the emergence of 
high-speed digital computers and automated 
modeling software. The realistic behavior of 
even extremely complicated systems can now be 
effectively modeled using finite element analyses 
and/or complex system simulations embodying 
relevant characteristics such as nonlinearities, 
coupling, time-varying parameters, large 
numbers of inputs and outputs, time delays, and 
input-derivative terms. Most models of complex 
systems are inherently high-order, even when the 
other complexities identified above are eliminated 
via linearization, estimation, or prudent 
disregarding of the offending terms. 
Consequently, the control system designer must 
often cope with unwieldy, high-order models 
which lead to the design of inordinately high- 
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order control systems. Although, in theory, the 
controller developed provides the desired 
performance, implementation of the full-order 
controller is sometimes not practical. 

In recent years, a considerable amount of 
research effort has been directed towards the 
problem of developing reliable order-reduction 
techniques for both mathematical system models 
and their associated controllers. Reduction in 
model order generally leads to a simpler, lower- 
order controller design, thereby increasing the 
feasibility and practicality of implementation and 
potentially reducing cost. The subsequent 
reduction in the degrees-of-freedom associated 
with the reduced-order model aids in the designer's 
visualization and understanding of the system and 
may lead to further insight into a control system 
design. Typically, the complexity of a model is 
reduced by first eliminating or simplifying 
nonlinearities, coupling effects, time-varying 
parameters, and other complicating factors. The 
model is then further simplified by reducing the 
order of the less complex model via frequency 
domain or state-space techniques. Control 
systems designed in this manner can provide 
acceptable performance in principle. However, in 
some instances the controller may require further 
tuning to improve performance, or may 
necessitate the addition of filters to eliminate the 
effects of unmodeled dynamic modes that interact 
with the control system in an undesirable 
fashion. 

The technique presented and illustrated 
in this paper constitutes a new and innovative 
way to achieve effective reduction of model- 
complexity by active control. The examples 
presented are simple engineering examples, 
however they demonstrate the potential for 
application of this method to more complicated 
systems. The traditional approach of first 
reducing the model-complexity and then 
designing the control system provides no 
mechanism whereby the credibility of the 


reduced-complexity model (RCM) is guaranteed. 
The term "reduced-complexity'* is used in place of 
"reduced-order" because the Model Complexity 
Reduction by Active Control (MCRAC) method 
produces a controller that is inherently capable of 
eliminating or minimizing the undesirable effects 
of phenomena such as coupling and nonlinear 
terms, as will be demonstrated via example 
problems. The method links the model- 
complexity reduction with the controller design 
method such that upon implementation of the 
controller, the credibility of the assumed RCM is 
automatically enforced by the controller actions, 
thereby achieving a new level of effectiveness in 
applying the classical linearization techniques. 

The proposed method employs 
"disturbance observer" theory to generate real- 
time estimates of the time-domain effects of the 
unwanted terms associated with model 
complexity. The total control effort is divided in 
such a way that part of the controller action 
annihilates, or minimizes, the disturbing effects 
of the model-complexity terms on the time 
response of the system. The other part of the 
control effort provides the desired closed-loop 
system performance. The control system 
designer is thus free to tailor the desired RCM 
and develop a control strategy such that the 
closed-loop RCM meets the desired performance 
goals and rejects disturbances acting on the 
system. The following sections present details 
of the method followed by three illustrative 
examples. Included in the Summary and 
Conclusions section are plans to extend the 
method to more complicated problems. 

Model Complexity Reduction 

Model complexity reduction is more 
generally referred to as model order reduction and 
is considered to be the reduction of the dynamic 
order of the model. Given a system of order n, 
the goal of the model order reduction process is 
to derive a model of order k, where k<n, such 
that for a given input, the reduced-order model 
output(s) closely tracks that of the full order 
system model. An error norm is examined to 
determine the accuracy of the reduced-order model. 
A plethora of methods exists by which this may 
be accomplished using time- and frequency- 
domain techniques as indicated in the extensive 
bibliographies in References 1 and 2. Indeed this 
area of model simplification has garnered the 
majority of the attention of researchers in this 
area. Techniques including balanced realization, 
component mode synthesis, aggregation 
methods, perturbation methods, and continued 


fraction methods have been developed to 
accomplish the task. 

Less frequently addressed are other 
aspects of model complexity reduction. These 
include approximation procedures for dealing 
with nonlinearities, time-varying parameters, 
coupling, and the number of system inputs. 
System nonlinearities are typically addressed by 
conventional linearization about a nominal 
operating point (or points) or a nominal 
trajectory. Functional linearizations such as 
describing functions may also be employed. 
Time-varying parameters are typically eliminated 
by using averaging methods or frozen-time 
eigenvalue methods. Coupling effects are 
sometimes ignored based on the assumption that 
they manifest fast transients. Coupling is also 
eliminated from high-order models using 
transformations to alternate coordinate systems, 
such as modal coordinates. The transformation 
masks the existence of coupling and the 
transformed model is then more easily reduced for 
controller design. For systems with a "large" 
number of inputs, often the number of inputs is 
reduced for simplification of the design via a 
fixed coordination of multiple control inputs or 
by ignoring certain disturbance inputs. Multiple- 
input multiple-output techniques are currently 
enjoying a much wider application to this 
modeling aspect. 

The MCRAC design philosophy is that 
a linear controller for the original full-complexity 
system model, M, can be designed using a 
judiciously chosen reduced complexity model, m, 
such that the resulting controller, u(m), provides 
performance which meets a given set of closed- 
loop specifications when implemented with the 
original system model, M. Furthermore, the 
action of u(m) enforces the RCM selected and 
guarantees the desired performance of M. The 
reduced complexity model is of lower order than 
the original and is linear. The resulting control 
design is, in general, of lower order than one 
designed using the original model. 

Assume that M has the standard state 
space form given in equation (1) below, where 
the matrix A may be time-varying, nonlinear, 
highly coupled, or embody other complexities. 

x = F(x,t,u,w) 

x = Ax + Bu + Fw + h.o.t.'s (1) 

y = Cx 
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In the equation above, x is the state variable, u is 
the control variable, w is an external disturbance, 
y is the output, and "h.o.t.'s" stands for higher 
order terms which arise from classical 
linearization. A reduced complexity model may 
be developed from (1) in the following manner. 
A linear RCM system matrix A r and the state 

vector x r are selected by the designer to embody 
the primary dynamics to be controlled. If x s is 
defined as the state vector which embodies the 
secondary dynamics, then the full system model 
can be written as given in equation (2) which, for 
now, ignores the h.o.t.'s derived from 
linearization. 


X = Ar A rs l/Xr\ + B u ( m ) + p w 

La„a,J (2) 

y = C r x T 

The reduced complexity model is now given by 


x r = A r x r + ArsXs + B r u(m) +Fw 


y = CrXr 


( 3 ) 


Control (DAC) methods. [3] Consider the 

external disturbance term and recall that for 
cancellation 


-B r iid(m) = Fw 4 (6) 

The disturbance w is modeled in state-space form 
as 


z = Dz + a (7) 

w = Hz. 

In equation (7), z is the state of the disturbance w 
and D is the system matrix governing the 
dynamics of z. The matrix D is developed 
directly from the time-domain or differential 
equations defining the waveform structure of w. 
The a are Dirac impulses which arrive at random 
intervals with random intensities. The derivation 
of the matrix D is explained thoroughly in the 
reference. From equation (6), disturbance 
cancellation occurs when 

Ud(m) = -B^FHz (8) 


which requires that 

with the dynamics of x s determined by equation 

(4) below. rank[B] = rankfBDFH], (9) 


X s = AsrX r + A,X, (4) 

Before explaining how the effects of the 
secondary dynamics are handled in the design, the 
structure of the controller u(m), given in equation 
(5) below, is examined. 

u(m) = Up(m) + Ud(m) +Uc(m) (5) 

where 

Up(m) = control providing desired response 
Ud(m) s control for external disturbances 
Uc(m) h control for system complexities. 


If the disturbance dynamics described in equation 
(4) manifest a waveform structure, which is 
typically the case for physical systems, then 
u c (m) can be designed using the same principles 
as described above for the external disturbances. 
A disturbance state model is developed based on 
the waveform structure of x s leading to the 
development of a state disturbance model as in 
equation (10). 

it = D,z, + a 

x, = Hz, (10) 

Therefore, for cancellation, 


Note that the total control effort u(m) is divided 
into three parts, the primary part u p (m) that 
provides the desired closed-loop system response, 
the disturbance accommodation part, utf(m) that 
is designed to cancel or minimize external 
disturbances, and the RCM part Uc(m) that is 
designed to cancel or minimize the effects of the 
unmodeled modes, nonlinearities (or h.o.t.'s 
resulting form linearization), and other 
complexities on the behavior of the reduced- 
complexity system. Those effects are viewed as 
time- varying disturbances, and thus u c (m) can be 
designed using Disturbance Accommodating 


Uc(m) = -Br A»Hz, ( 1 1 ) 

with the rank condition as indicated in equation 
(9) satisfied. When full cancellation occurs, the 
dynamic order of the system is reduced. The 
control Uc(m) must eliminate the effects of the 
coupling term, or nonlinear terms (demonstrated 
in the third example in the following section), 
without destabilizing the system. 

The total control action u(m) can be divided into 
additional parts, as required, to control any other 
disturbances that may act on the plant. The 
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up(m) part of the controller is designed using the 
engineer's method of choice so long as the 
controller meets the performance requirements. 
A composite observer is designed using standard 
observer theory to estimate the system states 
required for controller implementation as well as 
the effects of the unmodeled dynamics and 
disturbances. The observer order generally 
determines the order of the controller, therefore, 
care must be taken when selecting the disturbance 
model for the full-order model complexities. The 
choice of the disturbance model will be discussed 
in the following section. The outputs of the 
observer are used with the appropriate gains, 
calculated using the RCM, to implement u(m). 
The method will be demonstrated and elaborated 
upon using three simple examples in the 
following section. 

Illustrative Examples 

Three design examples are presented here 
to demonstrate the concept of the MCRAC 
method for controller design. It is emphasized 
that the results shown here are preliminary, but 
show a high degree of promise for the method. 
All examples are single-input, single-output 
systems with the assumption of co-location of 
the sensor and actuator. It is further assumed that 
the complexities removed from the full-order 
model may be modeled as an input disturbance. 
The first example comprises an eight-order 
spring-mass system. The second example is a 
highly-coupled sixth order system which 
simulates a flexible beam. The third example is 
a nonlinear second-order system. 

Shown in Figure 1 is a diagram 
representing an eighth-order dynamic spring-mass 
system, which is the focus of the first example. 
The values of the masses and spring constants 
used in this example are given below the diagram 
in Table 1. The control objective is to move mi 
to a given set point position and hold it there, 
regardless of the motion of the other three 
masses. Note that 



Figure 1. Spring-mass system for Example 1. 


m] = 10 

in 

II 

M 

m 2 = 5 

k 2 = 10 [ 

m 3 = 2.5 

II 

m 4 = 5 

M = 2.5 | 


Table 1. Parameter values used in Example 1. 


there is no damping in the system so that the 
masses m 2 , m 3 , and m 4 will continually 
oscillate. It is assumed that the position of mi 
can be sensed and that the control action is 
applied at mi. The equations of motion were 
derived and are given below in equations ( 12 ). 


y 1 _ J kl + k2 


yi + — Y2 + 

mi * 


m l 



y4 = ^-y3 - ^-y4 

ITI4 ITI4 


( 12 ) 


The eight-order model was cast in state space 
form for complexity reduction and control design. 
Because the goal of the closed-loop controller is 
to move mi to a particular position and hold it 
there regardless of the motion of the other 
masses, the reduced-complexity second order 
RCM, m, given in equations (13) below was 
selected. In this equation, y r represents the 
variable yi, therefore, the primary states for the 
system are those associated with yi. The 
motion of the other masses and the resulting 
forces imparted on mi are secondary and are 
modeled as a disturbance input w s$ which equates 
to x s in equation ( 10 ). 


y r + +-* 2. y r = -l-u(m) + — ^-w s 

mi mi mi 

(13) 

The control input u(m) is divided as in equation 
( 5 ) except that there is no ud(m) term as no 
external disturbances are considered. The primary 
control up(m) is designed such that the closed- 
loop response of m manifests the integral time 
absolute error (i.t.a.e.) response characteristics. 
The ideal design model is given in equation (14) 
below. 


y m + 1.4o) 0 ym + co 0 ym= y sp 


(14) 


4 



To design up(m), the standard set point error 
dynamics are examined for the ideal model and 
the RCM. When the error dynamics are equated, 
the gains for Up(m) are generated. The resulting 
controller providing i.t.a.e. response 
characteristics for the RCM is given below (for 
the parameter values given in Table 1). 


Up(m) = 10(g)o -1.5)(yt P -yi) - 14ooyi + 15y sp 

(15) 


For (00=10, 


Up(m) = 1000y sp - 985yi - 140yi (16) 

The character of the disturbance acting 
on mi is of waveform structure because of the 
sinusoidal motions exhibited by spring-mass 
systems. Therefore, the principles discussed in 
the previous section may be applied to develop 
the model of the disturbance for incorporation 
into the observer design. The motion of the 
system is sinusoidal, and certainly the 
disturbance may be modeled as a linear 
combination of sinusoidal functions of proper 
frequency. However, to be more general and to 
allow for unforeseen parameter variations, the 
spline disturbance model given in equation (17) 
is selected. The variable t is time and the q's are 
arbitrary constants. 

w s = c 1 t 2 + c 2 t + c 3 ( 17 ) 


The composite dynamics state model is given in 
equations (20). 



where 


A r - 


A ri H = 


0 

1 


(ki +k2\ 

0 

B r = 

m l / 


> 

0 0 0 

D, = 

010' 

^2. o 0 

001 

mi J 

J 

LoooJ 



c =[10000] 


The observer design equation is given in (21), 
where Xoi are the observer poles, placed for this 
design at X 0 i= -20, and Kq are the observer 
gains to be calculated. 

de{(A + KoC) - Xoil] = 0 (21) 


The resulting "disturbance" state model is given 
in (18) where D s and H are defined shortly. 

z, = D s z, + a 

w, = Hz. (18) 


For cancellation of the effects of the model 
complexities, the rank condition given in 
equation (11) must be satisfied, which is the case 
for this example. Therefore, the expression for 
uc(m) is derived immediately in equation (19). 


Uc(m) = -Br A ra Hz, 


Uc(m) 


“ 0 mi J" 


0 0 


Z, = -k2Zl 


^2. 0 
Lmi 


(19) 


An observer is designed using the 
composite dynamics involving m and the 
disturbance model for controller implementation. 


The design procedure for such observers is 
detailed in the DAC reference cited above. The 
resulting estimator, and thus controller, is 5th 
order. The numerical values of the gains for the 
state and disturbance estimates, indicated by the 
carets, for the closed-loop system are given in 
equation (22). Note that the disturbance 
derivatives are not required in the controller. 

u(m) - lOOOyip -985xi - 140x2 - 10zi ( 2 2) 

The plot in Figure 2 is a plot of the 
response of the full-order system model when 
controlled by u(m) as designed above. In the 
figure, the system output is the solid line, the set 
point is the dashed line, and the estimate of the 
output is another dashed line which tracks the 
true output exactly. Figure 3 is a plot of the 
disturbance (solid line), its estimate (tracks 
exactly the disturbance signal), and the control 
signal (dashed line). The control signal has been 
scaled in the plot to that it is evident that the 
controller is tracking the disturbance dynamics, 
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caused by the oscillations of the masses m2, 
m3, and 1114 and providing a signal which 
eliminates the effects of these oscillations on the 
position of mi. To test the robustness of the 
design to full-order model parameter variations, a 
series of simulations was executed while varying 
the mi’s and kfs for the system. The plot in 
Figure 4 shows that the system continues to 
track the set point input even with the system 
matrix parameters varied by 100%. No change to 
the control gains or the observer were made to 
achieve this result. 



time 


Figure 2. Set Point response with MCRAC 
controller. 



Figure 3. Control Signal and Disturbance. 

The second example, diagramed in 
Figure 5, focused on a simple model of a 
vertically suspended flexible beam. For 
simplicity, the beam is modeled as three rigid 
beams with torsional springs at the joints. 
Small angles are assumed in deriving the 
equations of motion for the sixth order system. 
The controller goal is to maintain 63 at zero 
degrees regardless of the motion of the other two 
rigid members. The system is highly coupled. 
The equations of motion for the system were 
derived using the Lagrange method and are 
presented in equations (23), with the definition of 


terms as indicated in figure 5 and the rj equal to 
the midpoint of each beam as measured from each 
joint. 



time 

Figure 4. Response with full-order model 
parameters varied by 100%. 



Figure 5. Rigid link "beam" model. 

©i(Ii+miri+(m2+m3)li) + 02(ni2lir2+ni3lil2) 

+ ©3m3hr3 + (mign+m2gli+m3gli)0t 
+ Ki0i +K 2 0i - K 2 02 = 0 

02(l2Hn2T2+ni3l2) + ©i(m2lir 2 +m3lil 2 ) 

+ ©3m3l2T3+ (m2gr 2 +m 2 gli+m3gl2)02 
+ K202 - K201 + K302 - K303 = 0 

03(l3+ni3r|) + 0im 3 lir3 
+ 02m3l2T3 + m 2 gr303 

-K302 + K303 = O (23) 

After obtaining a general expression for the 
mass and stiffness matrix, numerical values were 
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substituted and the full-order state space model M 
derived. For the results shown below, mj=12, 
Ij=l, kplO, and the length of each beam was set 
to 1. The RCM selected for this example was a 
simple second-order system involving the states 
for 03. The model took the same form as that for 
Example 1 . The same disturbance model was 
selected and the observer designed in a similar 
fashion to that in the previous example. The 
controller u(m), designed using m, was then 
implemented with the full-order model. In the 
simulation, the controller is activated at t=2 
seconds. The plot shown in Figure 6 is the ®3 
closed-loop response and the estimate, which 
tracks exactly. Figure 7 is a plot of the motion 
of 9i and 02 . 



Figure 6. Closed-loop response for Example 2. 


Example 2 Motion of Tbetol and Tbeta2 




time 


Figure 7. Motion of 01 (solid line) and 02 
(dashed line). 

Shown in Figure 8 is the control signal. 
Note that the controller produces a signal which 
cancels the effect of the motion of the upper two 
beam elements on the third. The controller is 
designed to accommodate all coupling effects and 
higher order terms as described in the previous 


section. These are estimated and eliminated 
through the MCRAC action. 



Figure 8. Control signal for Example 2. 


The third and final example involves the 
nonlinear second-order system given in equation 
(12). 

y + (l + y 2 )sin coy = 0 (12) 


The nonlinearity is the product of a sine function 
and the square of the derivative of the output y. 
The RCM chosen for this example was a simple 
double-integrator. An i.t.a.e. controller was 
designed for the RCM. The disturbance model 
selected to eliminate the effects of the nonlinear 
term is again a second-order spline implying a 
third-order disturbance model. The resulting 
controller is fifth order. The curves plotted in 
Figure 9 below indicate the system closed -loop 
response with u(m) in the loop with the full- 
complexity model. The initial conditions on the 

system are (y»y) = (1,0), the extremes on co are 
-100 < to < 100, and selected curves are plotted. 
The controller action is initiated at t=2 seconds. 
Similar results are obtained for varying the initial 
conditions on y and its derivative. 

A controller was designed using a 
linearized plant model for equation (12) and a 
similar set of simulation runs executed. The 
domain in which the controller, designed using 
the linearized model, produced acceptable results 
was -5 < to < 5.484. The controller designed 
using MCRAC ideas provides a substantially 
improved operating range over parameter 
variations. 
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Figure 9. Closed-loop response of nonlinear 
plant as to is varied from -100 to 100. 


2. Enns, D., "Model Reduction for Control 
System Design", Doctoral Thesis, Stanford 
University, June, 1984. 

3. Johnson, C. D., "Theory of Disturbance 
Accommodating Controllers", Chapter 7 in 
Advances in Control and Dynamics Sy stems, 
Vol. 12, (C. T. Leondes, ed.), Academic Press, 
New York, 1976. 


Summary and Conclusions 

The Model Complexity-Reduction by 
Active Control methodology has been introduced 
in this paper. The method shows potential for 
application to systems in which the assumptions 
described in the previous examples hold. That is, 
single-input, single-output systems with a co- 
located sensor and actuator pair. It is also 
assumed that the complexities removed from the 
full-order model may be modeled as an input 
disturbance. The controllers derived are of lower 
order than those derived using the full-complexity 
models for the two dynamic system examples and 
much higher order for the nonlinear system. 
However, the improvement in performance of the 
higher order controller over the second order 
controller designed using the linearized model for 
example 3 is significant. In all example cases, 
parameter variations of significant magnitude 
may be imparted on the full-order system matrix 
with virtually no effect on system performance. 

Future research will be directed toward 
further development of the theory and to the 
application of MCRAC principles to more 
complex systems for which the assumptions 
listed above do not necessarily hold. Based on 
the results presented in this paper, the method 
holds promise for application to the class of 
problems discussed. 
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